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ABSTRACT 

The spectrum observed by Fermi-LAT from the Galactic centre region shows a gamma 
ray feature near 130 GeV, that in some analyses appears as a possible line. We discuss 
the possibility that this gamma ray feature has a cosmic ray origin. We argue that 
the cosmic ray electrons steepen near 1 TeV from to about and are all 

secondary derived from the knee- feature of normal cosmic rays. We argue that the 
observed feature at ^ 130 GeV could essentially be a noise feature on top of a sharp 
turn-off in the 7 ray spectrum at ~ 130 GeV. This match suggests that the knee of 
normal cosmic rays is the same everywhere in the Galaxy. We suggest that it follows 
that all supernovae contributing give the same cosmic ray spectrum, with the knee 
feature given by common stellar properties; in fact, this is consistent with the super- 
nova theory proposed by Bisnovatyi-Kogan (1970), a magneto-rotational mechanism: 
Massive stars converge to common properties in terms of rotation and magnetic fields 
just before they explode. 

Key words: starsiwinds, outflows; gamma-rays:galaxies 



1 INTRODUCTION 

Fermi-LAT has revealed many interesting phenomena in the 
gamma ray sky. The Fermi bubbles have been observed ex- 
tending up to 10 kpc upwards and below from the Galac- 
tic Center (GC) (Su et al. 2010) in the 1-100 GeV energy 
range. More recently a Fermi group has reported the de- 
tection of GeV line emission from the Galactic centre re- 
gion (Ackermann et al. 2012). The diffuse gamma ray data 
taken by Fermi-LAT show a weak line or lines from 1.5° 
west of the Galactic Center (GC) (Su & Finkbeiner 2012). 
The observational features of the line like structure are not 
yet clear. More sensitive detectors are needed to confirm 
and resolve this structure (Li & Yuan 2012) and then it 
would be possible to know whether these are lines and dark 
matter signatures or some astrophysical broadband effect. 
Su & Finkbeiner (2012) have interpreted the observational 
data as a single spectral line at 127 ± 2 GeV or a pair 
of lines at 110.8 ± 4.4 GeV and 128.8 ± 2.7 GeV giving a 
slightly better fit. The analysis of the Fermi-LAT data by 
Cohen et al. (2012) prefers a single line either at 130 GeV or 



145 GeV. The line emission has been explained earlier with 
dark matter annihilations (Bringmann et al. 2012, Weniger 
2012, Tempel et al. 2012, Kyae & Park 2012, Lee et al. 
2012, Rajaraman et al. 2012, Buckley & Hooper 2012, Su 
& Finkbeiner 2012, Hea & Kim 2012, Yang et al. 2012). 

Cosmic rays may also explain the gamma ray feature at 
130 GeV. It has been suggested that a break in the gamma 
ray spectrum of the Fermi Bubbles may appear as an excess 
at 130 GeV (Profumo & Linden 2012). 

The interstellar total radiation field has been calculated 
by Moskalenko et al. (2006) as a function of the Galactocen- 
tric radius, taking into account the emission by stars of dif- 
ferent types and the scattering, absorption and re-emission 
of starlight by dust grains. It has distinct peaks and the 
photon density is maximum near 50-100/im and in the next 
peak near 1000/im. The intensity of the radiation field near 
the GC region is expected to be much higher than its locally 
observed value. In the present work we show that the 130 
GeV feature could be due to the Inverse Compton scattering 
of the infrared photons of wavelength 50/xm with the cosmic 
ray electrons of energy 1 TeV where the electron spectrum 
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steepens from E^'^ to about E^*"^ . We then show how this 
fits into a broader interpretation of recent cosmic ray data, 
and conclude that the underlying cosmic ray electron spec- 
trum near TeV energies very likely is the same in the GC 
region and near us. This strongly suggests that these elec- 
trons are secondary from spallation of cosmic ray nuclei near 
the knee energies, of about a few PeV; the cosmic ray nuclei 
spectrum then also has to be essentially the same at high 
energy in the Galactic centre region and near us, and quite 
specifically show the knee feature at the same energy. 



2 STAR FORMATION AT GALACTIC 
CENTRE AND FERMI BUBBLES 

Multi-wavelength observations have revealed surprising phe- 
nomena near the core of our Galaxy. The Fermi-LAT gamma 
ray detector has observed two large bubbles symmetrically 
located below and above the centre of our Galaxy. The bub- 
bles extend up to 10 kpc above and below the Galactic plane 
with a width of 40°. 

More interestingly, the WMAP (Jarosik et al. 2011 ) 
and Planck (Ade et al. 2012) haze is located within the 
northern bubble sharing the same edges and the ROSAT 
soft X-ray maps (Snowden et al. 1997) are in the circumfer- 
ence of the bubbles. These observational correlations suggest 
a common origin of these emissions. The gamma ray energy 
flux in the energy range of 1-100 GeV from Fermi Bubbles 
is E^dN{E^)/dE^ = 3 x 10"^ GeV cm"^ s"^ sr"^ and the 
solid angle subtended is 0.808 sr (Su et al. 2010). 

This gamma ray emission has been explained as a re- 
sult of cosmic ray protons interacting with protons in the 
medium (Crocker & Aharonian 2011) which leads to the pro- 
duction of charged and neutral pions. In this scenario if the 
cosmic ray protons are trapped for timescales of lO^^years, 
their energy partially going to the high energy gamma rays 
produced from decaying neutral pions can explain the ob- 
servations by Fermi-LAT. 

A more detailed discussion on mass and energy flows 
through the Galactic centre and into the Fermi Bubbles is 
given in Crocker (2012). Within a single-zone model it is 
shown that the star formation in the Galactic centre can 
sustain the energetics of the system and it explains the ob- 
served non-thermal emissions. 

Cheng et al. (2011) have interpreted the observational 
features near the core of our Galaxy as a result of periodic 
star capture processes by the galactic supermassive black 
hole SgrA* . Hot plasma wind is injected into the halo which 
subsequently heats up the halo gas. The hot gas produces 
thermal X-rays. The periodic injection of hot plasma also 
produces shocks and accelerate the electrons to TeV energy. 
Their synchrotron and Inverse Compton emission may ex- 
plain the observed radio and gamma ray emission from the 
central part of our Galaxy. 

The bubbles are also explained as evidence of possible 
AGN jet activity in our Galaxy (Guo et al. 2011), last quasar 
outburst (Zubovas et al. 2011) and coUimation of a wide an- 
gle outflow from SgrA* (Zubovas et al. 2012) by the Central 
Molecular Zone. 

The bubbles could be the result from second order Fermi 
acceleration of electrons by plasma wave turbulence through 
out the entire volume of the bubbles (Mertsch & Sarkar 



2011). The authors of this paper have explained the ob- 
served gamma ray emission with Inverse Compton losses of 
the hard spectrum of electrons resulting from second order 
Fermi acceleration. 

The magnetic field structure in the northern Fermi Bub- 
ble has been explored with polarized microwave emission 
(Jones et al. 2012). Their study reveals that the magnetic 
field lines in the northern Bubble's eastern wall and Galac- 
tic centre Spur are almost perpendicular to their extensions 
above the Galactic plane. 

We consider the case of a bubble being produced by the 
action of star formation in the Galactic centre region. Recent 
observations deduce a star formation rate of ~ 0.04-0.08 Mq 
yr"\ lasting about ~ 1-10 Myr (Yousef-Zadeh et al 2009; 
Immer et al 2012). The combined rate of momentum depo- 
sition from supernovae and stellar winds is then estimated 
to be F ~ 5 X 10^^{SFR/1 Mg/yr) ~ 2.5 x 10^^ dyne for a 
SFR ~ 0.05 Mq yr"^ (Leitherer et al 1999). There is also 
likely to be a component of pressure from cosmic rays driv- 
ing this wind (Everett et al 2008; Uhlig et al 2012). If we 
consider a bubble blown by the combined momentum injec- 
tion traversing a distance of L ~ 10 kpc in 15 Myr then the 
total energy deposited is~_FxL~7.7x 10^^ erg. A fraction 
~ 15% of this energy going to cosmic rays would produce a 
cosmic ray energy budget of ~ 10^* erg, consistent with the 
estimate of Su et al. (2010). The implied energy density in 
a bubble of radius ~ 5 kpc is ~ 1.4 x lO"^'^ erg cm~^, and 
the equipartition magnetic field is estimated to be ~ 4 fiG, 
consistent with the recent measurement of magnetic field in 
the bubble by Carretti et al. (2013). 



3 THE COSMIC RAY SPECTRUM IN THE 
GALACTIC CENTER REGION 

The spectrum of cosmic rays is determined by their injec- 
tion spectrum, their transport processes, and their losses. In 
a series of papers from 1993 (Biermann 1993; Biermann & 
Cassinelli 1993; Biermann & Strom 1993; Stanev et al. 1993) 
it has been proposed, that the injection of Galactic cosmic 
rays can be divided into two main source classes, one due 
to explosions into the interstellar medium, and one due to 
explosions into the predecessor stellar wind. For the ISM- 
SN-CRs the prediction had been a spectrum of £;-2-75±o.04 
after taking into account transport processes (with a Kol- 
mogorov spectrum), £-2-67-0.04 wind-SN-CRs below the 
knee, and the spectrum of £;-3 07-o.i4 g^j^Q^g ^j^g knee. The 
most recent fits to data give spectral index -3.08 beyond 
the knee (Biermann et al. 2003, Abbasi et al. 2012a, 2012b, 
Biermann & de Souza 2012). 

The knee energy scales with the charge of the nucleus, 
and is about 600 Z TeV, with a large uncertainty in its orig- 
inal prediction. The wind- Supernovae also have a polar cap 
component with an injection spectrum of E~^ and a sharp 
cutoff at the knee-energy due to a spatial limit in the accel- 
eration region. The final cutoff for ISM-SN-CRs had been 
predicted to be around 100 Z TeV, and for wind-SN-CRs at 
about 100 Z PeV; these predictions had been summarized in 
Biermann (1994), and Biermann et al. (2001). Since 2009 a 
series of tests have been possible that checked on the predic- 
tions, the CR-electron and CR-positron spectra (Biermann 
et al. 2009), the WMAP haze and 511 keV emission line 
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Figure 1. Observed electron spectrum near the earth propor- 
tional to and below and above the break at 1 ToV 
respectively. 



(Biermann et al. 2010a), the spectra of the various elements 
constituting the ISM and the wind components (Biermann 
et al. 2010b), and the different spallation histories for ISM- 
SN-CRs and wind-SN-CRs (Nath et al. 2012), with a low in- 
teraction column for ISM-SN-CRs and a high column for the 
wind-SN-CRs. The spectrum of the cosmic rays beyond the 
knee was tested in Biermann & de Souza (2012). These pre- 
dictions have also been used in explanations of the ultrahigh 
energy cosmic rays (Gopal-Krishna et al. 2010, Biermann & 
de Souza 2012), invoking a shift of the normal cosmic rays 
by a one-cycle relativistic shock following Gallant & Achter- 
berg (1999), and Achterberg et al. (2001) in the neighboring 
radio galaxy Cen A. The Fermi data of the Galactic Center 
region provide another test for these proposals. 

Here we address the possibility that there is a sharp 
spectral feature at 130 GeV in the region of the Galactic 
Center. For this we need to summarize the properties of 
the relevant cosmic ray electrons. The polar cap component 
of E~'^^^ in observer frame has a sharp cut-off at the knee 
energy, while the An component goes from E~^^^ to about 
E~^'^ at exactly the same energy. Pitch angle scattering and 
smoothing in the spallation production of the secondary lep- 
tons will modify this original spectral shape. After propa- 
gation the polar cap component of electrons from wind-SNe 
steepens to E~^~^^^ due to diffusion losses and at higher en- 
ergy E~^~^ due to synchrotron and Inverse Compton losses. 

Interestingly, a break near 1 TeV in the electron spec- 
trum has been observed near Earth (Biermann et al. 2009) 
and it is expected to be present also near the Galactic centre 
region. 

The knee near 4 PeV in the high energy cosmic ray 
spectrum can be correlated with the break in the electron 
spectrum at 1 TeV. Near the knee the cosmic rays are mostly 
He and heavier nuclei. The Lorentz factor of the heavy nuclei 
at the knee 'yknee ~ Z/A'jp, where 7p is the Lorentz factor of 
the protons. Z/A is nearly 1/2 for all heavy nuclei. We there- 
fore find ■jknee = 2 X 10®, which is also the Lorentz factor of 
the 1 TeV cosmic ray electrons at the break. The secondary 
electrons are perhaps produced in nuclear de-excitation fol- 
lowing energetic particle interactions (Ramaty et al. 1979) 



and as a result their Lorentz factors are approximately the 
same as that of the parent nuclei. 

The cosmic ray electron spectrum observed near Earth 
is shown in Figure 1: in the energy range of 100 GeV to 
4 TeV. The solid line shows the electron spectrum propor- 
tional to E~^ and E~'^'^ below and above the break at 1 TeV 
respectively. We will show below that if the steepening from 
E~'^ to E~*'^ is at 1 TeV then the Inverse Compton scat- 
tering of the infrared background photons by the electrons 
at the break energy can explain the 130 GeV line feature in 
the 7 ray spectrum near GC. This requires the cosmic ray 
electron spectrum around TeV energies to be essentially the 
same in shape near the Galactic Center as near Earth. 



4 130 GEV FEATURE EMISSION NEAR GC 

The interstellar radiation (IR) field has been estimated 
(Moskalenko et al. 2006) at different distances from the GC. 
The intensity of the field is highest at GC and decreases with 
distance from GC. The polar cap component of the electron 
spectrum (in the loss dominant case) as observed near the 
Earth is 

dNe K 



23.55 X 10"S^^cm"^ 



d'ye 

We assume near the Galactic centre region the electron flux 
is 77 > 1 times the flux observed near the Earth. The 7 
ray flux produced by Inverse Compton scattering of infrared 
photons near GC region is calculated. 

Vol dEic 1 



E^/jfm = E 



dE^ 



'''in Die dVdE-ydt AD. 



We are considering a region with —15° ^ & ^ 15° and 
—20° ^ / ^ 20° near the Galactic centre, which is similar 
to the region 3 given in Fig.l. of Weniger 2012. The solid 
angle subtended by this region to the observer on the Earth 
is AQ = 2 X d(f> j^Qo_^c^o sinOdO = 0.36sr, where ((> — I 
and 9 = 90° — b. The distance to GC from the Earth is 
Dec = 7.5 kpc (Nishiyama et al. 2006, Eisenhauer et al. 
2005) which gives AtvDIc = W*^'^ cm^ 

The power emitted per unit volume, per unit energy and 
per unit time by the gamma rays produced in Inverse Comp- 
ton scattering is expressed as ^YdE^dt • ^^^^ electron spectral 
index p = 3 gives gamma ray spectral index (p— l)/2-f 1 = 2, 
''^^f "'^ E-y^- The expression for power emitted in Inverse 
Compton scattering by electrons is (Rybicki & Lightman 
1979) 



dEic 
dVdE-ydt 



= 23.55 X lO^SycvrroAp^B; 



ev{e)de (3) 



where w(e) is the low energy photon density per unit energy 
e. The value of the constant Ap for p = 3 is 16/9, the classi- 
cal electron radius rg = 2.817 x 10"^'^ cm and c is the speed 
of light. We have taken rj x 1/0^/(10®"*'® cm^) ~ 10, where 77 is 
taken to be of order unity in order to take into account the 
possible variation in CR flux in the galactic centre area, and 
also the uncertainty in our knowledge of the volume. After 
taking angular average the electron Lorentz factor (7^), tar- 
get photon energy (e) and scattered photon energy (E-y) are 
related as 4/37^6 = Ej. Thus the electrons of energy 1 TeV 
Inverse Compton scatter the IR photons of energy 0.024 eV 
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Figure 2. Gamma ray data with error bars from Rcg3 in Fig.l 
of Weniger 2012, solid line sliows spectral shape with E~'^ and 
g— 2.6 below and above 130 GeV respectively. 

which corresponds to wavelength 50/im and 130 GeV 7 rays 
are produced. 

The IR background has several peaks as shown in Fig- 
ure 1 of Moskalenko et al. (2006). The IR photon density 
at the peak near 0.024-0.012 eV is higher than the densities 
at the higher energy peaks near 1.2 eV and 0.12 eV. It is 
comparable to the photon density at the next lower energy 
peak at 0.0012 eV. The IR photons of energy 0.024 eV are 
mostly contributing to the gamma ray spectrum near 130 
GeV. The lower energy IR photons produce gamma rays of 
the same energy by Inverse Compton scattering off higher 
energy electrons. As the electron spectrum steepens beyond 
1 TeV the contribution to the gamma ray flux at 130 GeV 
from the lower energy IR photons near 0.0012 eV is not sig- 
nificant. We have parametrized the IR peak near 0.024 eV as 
a delta function at that energy. We have used IR photon den- 
sity 1/(0.024)2 eV"^ cm ^ in our calculations (Moskalenko 
et al. 2006). The flux of the gamma rays is calculated us- 
ing eqn((2]), eqn((3]) and compared with the observed flux in 
Figure 2 of this paper. 

The calculated gamma ray spectrum is proportional to 
E~^'^ above 130 GeV. The sharp decrease in photon flux 
near 130 GeV due to the change in the spectral index of the 
electron spectrum has perhaps been observed as a line-like 
feature in the Fermi LAT data. 



5 DISCUSSIONS 

The gamma ray emission observed by Fermi-LAT near the 
GC region has been explained as Inverse Compton emission 
of cosmic ray electrons. The IR background photons in this 
region are Inverse Compton scattered off the cosmic ray elec- 
trons. The Galactic wind driven by star formation in the GC 
region and the emitted cosmic rays supplies the energy of 



~ 10^ erg to the bubbles of gamma rays. The break in the 
bubble's spectrum due to the break in the spectrum of the 
cosmic ray electrons appears as a line near 130 GeV. 

The intensity of gamma rays from the bubbles does not 
show any significant variation with distance from GC. This 
may be understood as follows: a) The gamma ray intensity 
of Fermi Bubbles depends on the intensity of the background 
IR radiation field and the density of the cosmic ray electrons. 
The IR radiation is essentially constant initially up into the 
halo as long as the lateral scale of the emitting region is 
larger than the distance above the disk, b) The density of the 
cosmic ray electrons is inversely proportional to the cross- 
section of the fiow. If the initial flow is straight up, then 
the density changes only slowly and adiabatic losses may be 
weak. Moreover, re-acceleration by weak shocks (eqn. 2.45 
in Drury 1983) may keep the spectral shape of the cosmic 
ray electrons intact but sharpen the kink in the spectrum, 
c) What we observe is a line of sight integral, so even a 
decrease of the IR radiation field weaker or similar to 
would be compensated in the integral for a straight initial 
flow up/down. Taken together this allows to understand the 
constancy of the emission. 

Fermi-LAT has also observed emission of gamma rays 
in 1-100 GeV band in the star forming region of Cygnus 
X (Ackermann et al. 2011). The observed gamma ray spec- 
trum is proportional to E~^. The polar cap component of 
cosmic rays at the source is proportional to and it gives 
secondary pions in hadronic interactions. The neutral pions 
subsequently decay to gamma rays. Thus the gamma ray 
spectrum observed from Cygnus X might have originated in 
hadronic interactions of the polar cap component of cosmic 
rays. IC emission spectrum of secondary electrons in this 
case would be proportional to E^^^'^ and inconsistent with 
the observed spectrum E^^ . 

The direct observations of the cosmic ray electron spec- 
trum near the Earth give us the loss-dominant spectrum 
of the cosmic ray electron spectrum in our neighborhood, 
and the gamma-ray observations give us that spectrum in 
the Galactic Center region. It is very interesting that the 
derived spectral shape is the same. 

It is possible to relate the knee in the very energy cos- 
mic ray spectrum with the break in the cosmic ray electron 
spectrum. Cosmic ray interactions and subsequent nuclear 
de-excitation leads to the production of the secondary elec- 
trons, as the cosmic ray parent nuclei and secondary elec- 
trons have the same Lorentz factor. This implies that the 
cosmic ray knee is at the same energy in terms of E/Z ev- 
erywhere in the Galaxy as the break energy in the electron 
spectrum is same in different regions of the Galaxy, as pre- 
dicted from a theory, that attributes this kink to the original 
stars. This derives from the fact, that in different parts of the 
Galaxy very different stars contribute to the kink feature of 
cosmic rays, the knee, and in each location very many stars 
contribute. So to have a clear feature at all, and then to have 
that feature the same in different parts of the Galaxy puts 
very strong constraints on the exploding stars: they must all 
be asymptotically very similar at the point of explosion. 

We need to ask, whether any other mechanism could 
produce such a kink in the cosmic ray spectrum, and also 
give a break energy which is the same everywhere in the 
Galaxy: First we may consider the possibility that OB super- 
bubbles produce the cosmic ray component giving the knee; 
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in such a theory a curvature upwards and below the energy 
of the knee would be consistent with arguments on shock 
structure (e.g. Drury 2011); One might speculate that it 
could also produce a knee from a characteristic length scale 
of the super-bubble. However, it is hard to see how this crit- 
ical length scale can be the same everywhere in the Galaxy. 
Second, we can generalize this argument to the transport 
in the Galaxy, and the escape from the Galaxy. However, 
as is well known, this possibility would give a much larger 
cosmic ray anisotropy in arrival directions than what has 
been observed, and so is also wrought with difficulties (e.g. 
Biermann 1993, and the references mentioned there). 

So, it appears that the original exploding stars are re- 
quired to explain this kink, called the knee. This has actu- 
ally been predicted (Biermann 1993), and is based on the 
magneto-rotational mechanism for the explosions of very 
massive stars proposed by Bisnovatyi-Kogan (1970), and 
worked out in much more detail by Ardeljan et al. (2005), 
Moiseenko et al. (2006), Bisnovatyi-Kogan & Moiseenko 
(2007) and others. This mechanism connects the rotation 
of the core of the star with the magnetic fields, just as is 
required to explain a constant knee E/Z scale (Biermann 
1993). 

We can also make further tests of the proposal. If the 
segment of the cosmic ray electron spectrum is sec- 
ondary, then the upturn should be connected to the paral- 
lel upturn in the positrons (e.g. Biermann et al. 2009, now 
Adriani et al. 2011); in fact this is consistent. 

The cosmic ray electrons are also emitting synchrotron 
photons which have been observed in the Planck spectrum 
(Ade et al. 2012). One zone compatible with radio spectrum 

due to the electron spectrum _E^'' in the loss limit and 
another with a spectrum compatible with u^^^^ due to the 
electron spectrum E^'^^'^ in the diffusion limit are expected 
to be observed by Planck. A ring may be present at the 
transition region between the outer zone of high frequency 
radio spectrum proportional to and inner zone where it 
is proportional to ly"^^'^ . Perhaps due to ignoring the ring 
region the Planck data is at present consistent with a slightly 
steeper spectrum in the outer zone and a flatter spectrum 
in the inner zone. 

It is also expected that there should be many wind-SN- 
remnants with a spectrum of E~'^^^ at lower gamma ray 
energies, and at higher energies. 



6 CONCLUSIONS 

In this work we propose an explanation of the 130 GeV line- 
like feature appearing in the GC gamma ray emission in 
the Fermi data: The idea is that this is just a noise feature 
on top of a sharp kink in the underlying continuum. The 
continuum is explained as Inverse Compton emission of the 
FIR radiation field by secondary cosmic electrons/positrons 
produced in spallation of cosmic ray nuclei at the knee. 

We realize that this implies that the cosmic ray elec- 
tron spectrum has the same shape at high energy in the GC 
region as near to Earth in the Galaxy. This implies then in 
our concept that the cosmic ray particles themselves have 
the same knee energy through out the Galaxy. This in turn 
supports an origin of this feature in the original supernova 



explosions, which in turn supports the magneto-rotational 
mechanism proposed by Bisnovatyi-Kogan (1970). 

If all this is true, we have identified a common quan- 
titative feature of all very massive stars in the final stage 
just before they blow up, connecting rotation and magnetic 
fields. 
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